Numerical modelling of thermal effects in rats due to high-field magnetic resonance imaging (0.5-1 GHz)

A finite-difference time-domain (FDTD) thermal model has been developed to compute the temperature elevation in the Sprague Dawley rat due to electromagnetic energy deposition in high-field magnetic resonance imaging (MRI). The field strengths examined ranged from 11.75-23.5 T (corresponding to H-1 resonances of 0.5-1 GHz) and an N-stub birdcage resonator was used to both transmit radio-frequency energy and receive the MRI signals. With an in-plane resolution of 1.95 mm, the inhomogeneous rat phantom forms a segmented model of 12 different tissue types, each having its electrical and thermal parameters assigned. The steady-state temperature distribution was calculated using a Pennes 'bioheat' approach. The numerical algorithm used to calculate the induced temperature distribution has been successfully validated against analytical solutions in the form of simplified spherical models with electrical and thermal properties of rat muscle. As well as assisting with the design of MRI experiments and apparatus, the numerical procedures developed in this study could help in future research and design of tumour-treating hyperthermia applicators to be used on rats in vivo.

A distributed equivalent magnetic current based FDTD method for the calculation of E-fields induced by gradient coils

This paper evaluates a new, low-frequency finite-difference time-domain method applied to the problem of induced E-fields/eddy currents in the human body resulting from the pulsed magnetic field gradients in MRI. In this algorithm, a distributed equivalent magnetic current is proposed as the electromagnetic source and is obtained by quasistatic calculation of the empty coil's vector potential or measurements therein. This technique circumvents the discretization of complicated gradient coil geometries into a mesh of Yee cells, and thereby enables any type of gradient coil modelling or other complex low frequency sources. The proposed method has been verified against an example with an analytical solution. Results are presented showing the spatial distribution of gradient-induced electric fields in a multi-layered spherical phantom model and a complete body model. (C) 2004 Elsevier Inc. All rights reserved.

Electromagnetic fields inside a lossy, multilayered spherical head phantom excited by MRI coils: models and methods

The precise evaluation of electromagnetic field (EMF) distributions inside biological samples is becoming an increasingly important design requirement for high field MRI systems. In evaluating the induced fields caused by magnetic field gradients and RF transmitter coils, a multilayered dielectric spherical head model is proposed to provide a better understanding of electromagnetic interactions when compared to a traditional homogeneous head phantom. This paper presents Debye potential (DP) and Dyadic Green's function (DGF)-based solutions of the EMFs inside a head-sized, stratified sphere with similar radial conductivity and permittivity profiles as a human head. The DP approach is formulated for the symmetric case in which the source is a circular loop carrying a harmonic-formed current over a wide frequency range. The DGF method is developed for generic cases in which the source may be any kind of RF coil whose current distribution can be evaluated using the method of moments. The calculated EMFs can then be used to deduce MRI imaging parameters. The proposed methods, while not representing the full complexity of a head model, offer advantages in rapid prototyping as the computation times are much lower than a full finite difference time domain calculation using a complex head model. Test examples demonstrate the capability of the proposed models/methods. It is anticipated that this model will be of particular value for high field MRI applications, especially the rapid evaluation of RF resonator (surface and volume coils) and high performance gradient set designs.

Winding patterns for biplanar MRI shim coils with rectangular and circular target-field regions

A method is presented for calculating the winding patterns required to design independent zonal and tesseral biplanar shim coils for magnetic resonance imaging. Streamline, target-field, Fourier integral and Fourier series methods are utilized. For both Fourier-based methods, the desired target field is specified on the surface of the conducting plates. For the Fourier series method it is possible to specify the target field at additional depths interior to the two conducting plates. The conducting plates are confined symmetrically in the xy plane with dimensions 2a x 2b, and are separated by 2d in the z direction. The specification of the target field is symmetric for the Fourier integral method, but can be over some asymmetric portion pa < x < qa and sb < y < tb of the coil dimensions (-1 < p < q < 1 and -1 < s < t < 1) for the Fourier series method. Arbitrary functions are used in the outer sections to ensure continuity of the magnetic field across the entire coil face. For the Fourier series case, the entire field is periodically extended as double half-range sine or cosine series. The resultant Fourier coefficients are substituted into the Fourier series and integral expressions for the internal and external magnetic fields, and stream functions on both the conducting surfaces. A contour plot of the stream function directly gives the required coil winding patterns. Spherical harmonic analysis of field calculations from a ZX shim coil indicates that example designs and theory are well matched.

An FDTD model for calculation of gradient-induced eddy currents in MRI system

In most magnetic resonance imaging (MRI) systems, pulsed magnetic gradient fields induce eddy currents in the conducting structures of the superconducting magnet. The eddy currents induced in structures within the cryostat are particularly problematic as they are characterized by long time constants by virtue of the low resistivity of the conductors. This paper presents a three-dimensional (3-D) finite-difference time-domain (FDTD) scheme in cylindrical coordinates for eddy-current calculation in conductors. This model is intended to be part of a complete FDTD model of an MRI system including all RF and low-frequency field generating units and electrical models of the patient. The singularity apparent in the governing equations is removed by using a series expansion method and the conductor-air boundary condition is handled using a variant of the surface impedance concept. The numerical difficulty due to the asymmetry of Maxwell equations for low-frequency eddy-current problems is circumvented by taking advantage of the known penetration behavior of the eddy-current fields. A perfectly matched layer absorbing boundary condition in 3-D cylindrical coordinates is also incorporated. The numerical method has been verified against analytical solutions for simple cases. Finally, the algorithm is illustrated by modeling a pulsed field gradient coil system within an MRI magnet system. The results demonstrate that the proposed FDTD scheme can be used to calculate large-scale eddy-current problems in materials with high conductivity at low frequencies.

Novel target-field method for designing shielded biplanar shim and gradient coils

A new method is presented here for the systematic design of biplanar shielded shim and gradient coils, for use in magnetic resonance imaging (MRI) and other applications. The desired target field interior to the coil is specified in advance, and a winding pattern is then designed to produce a field that matches the target as closely as possible. Both gradient and shim coils can be designed by this approach, and the target region can be located asymmetrically within the coil. The interior target field may be matched at two or more interior locations, to improve accuracy. When shields are present, the winding patterns are designed so that the fields exterior to the biplanar coil are made as small as possible. The method is illustrated here by the design of some transverse gradient and shim coils.

A diffusion and T2 relaxation MRI study of the ovine lumbar intervertebral disc under compression in vitro

The ovine lumbar intervertebral disc is a useful model for the human lumbar disc. We present preliminary estimates of diffusion coefficients and T-2 relaxation times in a pilot MRI study of the ovine lumbar intervertebral disc during uniaxial compression in vitro, and identify factors that hamper the ability to accurately monitor the temporal evolution of the effective diffusion tensor at high spatial resolution.

MR Microscopy and microspectroscopy of the intact kidney

Magnetic resonance (MR) techniques have become essential in clinical practice and a broad range of research areas. We begin with a review of the potential and limitations for resolution improvements by MR techniques. The kidney has distinct regional structural, functional and biochemical variability. The isolated perfused rat kidney (IPRK) retains renal function while eliminating movement and susceptibility boundaries which severely limit the potential of MR techniques. The IPRK, with a length of less than 20 mm in the longest axis, will be used to illustrate the potential resolution of different MR techniques and the different: biological information that can be obtained. (C) 2004 Wiley Periodicals, Inc.